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I. INTRODUCTION 



In previous work M, one of the present authors proposed a discrete singular convolution (DSC) algorithm for 
computer realization of singular convolutions involving singular kernels of delta type, Abel type and Hilbert type. 
One of illustrations for the algorithm was Shannon's sampling formulae Q which plays an important role in the 
approximation of the delta distribution and generalized derivatives ||. However, in practical computations, the 
truncation error of Shannon's sampling formulae is substantial . A regularization technique Q was used to construct 
a regularized Shannon's sampling formulae Jl]]^] , which was found to be extremely accurate and robust for resolving 
various challenging dynamical problems, such as the homoclinic orbit excitation of the Sine-Gordon equation (7), the 
Navier-Stokes flow in complex geometries p], shock capturing of the inviscid Burgers' equation |9], molecular quantum 
system described by the Schrodinger equation jlo| and nonlinear pattern formation of the Cahn-Hilliard equation . 

The objective of the present note is to provide a theoretical analysis for the previous excellent numerical results. 
Rigorous error estimations of the regularized Shannon's sampling formulae are given for their applications to inter- 
polations and derivatives of a function. 
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where superscript, (s), denotes the sth order derivative. 
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III. PROOF 



A. Separation of the error 



The error breaks naturally into a few components. Denote 
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Here, Ei(t) is regularization error. E 2 (t) and E 3 (t) are truncation errors. From Shannon's sampling theo 
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which can be differentiated term by term, the total error can be written as a sum of three components 
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The corresponding error norms satisfy 
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B. Estimation of E\(t) 

Let f(uj) be the Fourier transform of f(x), and f(u>) = J R f(x) e~xp(ixu))dx. Since 
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Since function / satisfies 
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it has a Fourier series expansion 
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then combining Eqs. (§), @, © and @, one has 
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For a; £ [— B, B], e(w) can be evaluated as 
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Moreover, for x > 0, the following inequality ]T2| is valid 
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Therefore, the estimation for e(w) is obtained as 
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It follows from Eqs. <M) and (E3J) that 
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From the Parseval identity, one has 
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C. Estimation of Eht(t) 



Differentiations can be written 
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where Hk(x) is the Hermite polynomial 
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Two simple lemmas are required. 

Lemma 1 (Abel's inequality) jL3|. For two sequences {a„}, {b n }, b\ > 62 > ■ • • > b n , a n , b n € -R, set 
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Lemma 2. As all the notations unchanged, set g(x) = f(x)Hk(^^)- If 5(2;) satisfies 

g'(x)<g(x)^-, 

whenever x > t + (Mi — 1)A, then {Fk(l)}ieAT is a decreasing sequence. 
The proof is obvious by taking the first order derivative. 
Let denote 
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Then from (28) and (37), and by using lemma 1 and lemma 2, one has 
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D. Estimation of E 3 (t) 

A result like lemma 2 is required. 

Lemma 3. Notations are the same as before. Denote g(x) = f(x)Hk{^^), if g(x) satisfies 

g'(x)>g(x)^-, 



whenever x < t — M2 A, then {Fk(l)}i^ is an increasing sequence. 

The proof is also direct. Therefore, by the same treatment as that in the previous subsection, we obtain 
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E. The end of the proof 

By combining Eqs. @, @, @ and (glj), one obtains Eq. (g). 



IV. DISCUSSION 



Remark 1. For s = 0, one has 
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This is a rigorous error statement for the formulae widely used in the aforementioned numerical computations. Roughly 
speaking, if exp(— ^-) = 10~ v , then rj = J" 10 , so the error is 
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where r 



One may choose r,B,A and M appropriately to attain desired accuracy. Assume all non-exponential 



quantities are combined to give unit, and M = Mi — 1 = M2, one has 



-(tt - BA) > v/?72 In 10 



and 



M 



— > yjr]2 In 10, 



(44) 



(45) 



where 77 is the desired order of accuracy. There are some general rules for attaining high accuracy. These are discussed 
from two different arguments. 

1.) For a given function f(x) with a known bandlimit £?, other parameters, A,r and M, are to be chosen appro- 
priately to achieve a desired accuracy order rj: 

V172 in 10 ^ -p Qr £ xe( j r ^ ^ ne ^ig^gj. the frequency bandlimit B is, the 



(i) From Eqs. (|4J) and ([45j) one has BA < it 
smaller A should be, which means the more grid points in the computational domain. When A varies from to -jy, r 

changes from ^ r)2 ln 10 to +00, therefore for sufficiently small A, r is near ^'' 2 ln 10 . 

(ii) No matter how many grid points are in the computational domain, r and M cannot be too small. Equations 



(||) and (|45|) indicate 
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and M > r\Jr\2 ln 10. If M and r are less than the minimal requirements, the 
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accuracy deteriorates quickly. On the other hand, if sufficiently large r and M are used, say, M = 30 and r = 3.5, 
high approximation accuracy can be achieved. 

2.) In practical computations, such as in solving a partial differential equation, the function / and its B are 
unknown. In this case, A is selected a priori. Then r and M are to be chosen properly for achieving a desired 
accuracy order 77: 

(i) For a given grid spacing, A, and accuracy requirement 77, r value determines frequency bandlimit B which can 
be reached. Then the set of functions / which are almost bandlimited to B can be accurately approximated (where 
'almost bandlimited to £?' means the function / is not necessarily bandlimited but its Fourier amplitude outside 
\B\ is much smaller than the given error 10^''). The choice of M should be consistent with r for a given accuracy 
requirement. In general, small r and M values lead to an accurate approximation for low frequency component of a 
function of interest. But the prediction of a high frequency component will not be accurate in such a case. 

(ii) For a given grid spacing A and r value, the larger M is, the higher bandlimit B can be reached. 

(iii) To improve computational efficiency with a given A, B shall be very close to ^. However, to maintain certain 
approximation accuracy, r has to be sufficiently large, which implies that M has to be very large too. This in turn 
results in low efficiency (It takes M — ► 00 to maintain the accuracy if one samples at the Nyquist rate) . 

(iv) If A, M and 77 are chosen, then r is also fixed. For example, to achieve the machine precision 10 _J? ~ 10~ 15 , 
Eq. (Ejj ) estimates r > 2.8. If this is achieved by using M = 33, then Eq. (|45| ) estimates r < 4. In fact, M ~ 30 and 
2.8 < r < 4.0 are the parameter regions found from an earlier numerical test |lj and were used in many applications 

Remark 2. A comparison between the truncation errors of Shannon's sampling formulae and the regularized 
Shannon's sampling formulae is in order. Reference 0] estimates that the expression 



has error of 



\T N (t)\ < ^VE 

TT 



■ f nt \ 

sm( A } 



^ (47) 



(TV A 2 —t 2 Y 

where t < NA, and E is the total 'energy' of the function given by 

This is not directly comparable with our error estimate because our sampling is centered around a point of interest, 
x. Let consider a truncation error of the form 

(E M f)(t) = m - £ /(^ y^ . (49) 



In Appendix A, it is shown that in a finite computational domain, the L 2 norm of (Emf)(t) has the order of 
\\f(t)\\L 2 (R)\J~jf%, which is much larger than the truncation error of the regularized Shannon's formulae. On the 

other hand, to achieve the same accuracy, the regularized formulae requires much fewer computational grids Jl|,^| . 

Remark 3. Discussions for the higher order derivatives can be presented in a similar manner as those of Remarks 1 
and 2. In fact, previous work of solving partial differential equations Jl|J^-pl]|. involved such derivatives, and results 
are consistent with the present theorem. Detailed comparison is omitted. 

Remark 4- In many practical applications, such as in solving partial differential equations, error estimations and 
discussions in other spaces are often required. Moreover, in real computations, the computational domain is always 
limited to a finite interval, such as [a, b]. Therefore, the norm ||/||x,2(^ in Eqs. ( |39| ) and ( |4l| ) are required to be 
changed into ||/|| L 2 (a,b)> wm ch can be evaluated by integrations along [a + M%A, b + Mi A] and [a — M2A, b — M2A] 
respectively. Therefore various L p (l < p < 2) error estimates of E(t) can be derived accordingly. If we know the size 
of LP norm (1 < p < 2) of the function of interest, then we can deduce from the theorem the critical values, r and 
M, to achieve desired accuracy. 
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APPENDIX A: TRUNCATION ERROR OF SHANNON'S SAMPLING FORMULAE 
Lemma 4. In the computational domain [a, b], the error expression 
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Proof. Let denote 
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By using Schwartz's inequality, one has 
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Similarly one obtains 
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By combining Eqs. (A3), (A4) and (A5), one finishes the proof. 



